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ABSTRACT.  The  need  to  accurately  model  multiscale  phenomena  is  ubiquitous  in 
eddy-current  nondestructive  evaluation.  By  using  volume-integral  equations,  we  are 
able  to  develop  a  very  simple  algorithm  for  accurately  computing  the  response  of  a 
very  small  anomaly  in  the  presence  of  a  much  larger  one.  We  validate  the  algorithm 
and  its  associated  code  in  VIC-3D©  through  benchmark  data  on  two  test  sets:  (1) 
a  notch  at  a  bolt  hole  with  an  upper  surface  coil,  and  (2)  a  notch  in  a  bolt  hole  with 
a  plate  surface  coil. 

Keywords:  Eddy-Current  Nondestructive  Evaluation,  Volume-Integral  Equations, 
Electromagnetic  Forward  and  Inverse  Problems,  Multiscale  Phenomena 
PACS:  02.50.-r,  81.70.-q 

INTRODUCTION 

The  use  of  integral  equations  and  anomalous  currents  allows  us  to  efficiently  re¬ 
move  ‘background  effects’  in  either  forward  or  inverse  modeling.  Consider  the  anoma¬ 
lous  region  within  a  background  host,  as  shown  in  Figure  1.  Let  <vb(r)  be  the  conduc¬ 
tivity  when  the  flaw  is  absent.  Outside  the  background  region,  <rb(r)  is  equal  to  the 
host  conductivity,  ah-  Inside  the  background  regions  it  varies  with  position,  r.  Let 
cr(r)  be  the  conductivity  when  the  flaw  is  present.  Outside  the  background  region,  <r(r) 
is  equal  to  the  host  conductivity.  Inside  the  background  region,  but  outside  the  flaw, 
<r(r)  =  <Tfe(r).  Inside  the  flaw,  <x(r)  is  not  equal  to  ab(r),  but  varies  with  position. 

First,  consider  the  unflawed  background  region.  The  anomalous  current  for  this 
problem  satisfies 

Jb(r)  =  (crb(r)  -  ah) (Ein(r)  +  E(r)[Jb])  ,  (1) 

where  the  functional  notation,  E(r)[Jb],  implies  an  integral  operator  on  Jb,  whose  kernel 
is  a  Green’s  function  [1], 

Next,  consider  the  flawed  background  region,  and  define  anomalous  currents 

J/(r)  =  (cr(r)  —  <Jb(r))(Ejn(r)  +  E(r)[J/])  (2) 

Jd(r)  =  {ah(r)-ah)(  Ein(r)  +  E^)^]  +  E(r)[Jd]) 

+  0(r)  -<76(r))E(r)[Jd]  .  (3) 
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FIGURE  1.  Defining  an  anomalous  region  within  a  background  host.  <75 (r)  —  ah  is  the 
anomalous  conductivity  associated  with  the  ‘large’  background,  and  <r(r)  —  c&(r)  is  the  anoma¬ 
lous  conductivity  associated  with  the  ‘small’  flaw. 


The  anomalous  current,  Ja  =  3d  +  J/,  satisfies  3a 
E(r)[JJ). 

The  change  in  impedance  due  to  the  flaw  is 


(cx(r)  -  ah)(Ein  +  E(r)  [J f]  + 


Ein  •  3 a  Ein  •  3b  Ejn  •  (J/  T  3d  J&) 

Ein  '  3 f  T  Ejn  '  Jint  j 

where  we  have  defined  Jm#  =  3d  —  J&. 

From  (l)-(3)  we  obtain  the  uncoupled  integral  equations 

Jft(r) 


Ezn(r) 

cr6(r)  -  ah 


Ein(r)  = 
Ee/(r)  = 


J/(r) 


afrj  - 


-  E(r)  [J6] 

-  E(r)  [Jf] 


-Er) 


Tint 


cr{ r)  -  ah 

where  the  effective  incident  field.  Eej(r),  is  given  by 

E«/(r)  =  [j  ]  +  g(r>  -  g»(r>  E  [JJ 

a(r)-ah  a(r)  -  ah 

The  algorithm  for  eliminating  the  background  is: 


1.  Solve  (5)  for  Jb  using  a  coarse  background  grid,  Gb 

2.  Solve  (6)  for  3f  using  a  fine  grid,  Gf,  covering  only  the  flaw 

3.  Solve  (7)  for  3int 


(4) 

(5) 

(6) 

(7) 

(8) 
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FIGURE  2.  Benchmark  Problem  1:  Notch  in  a  bolt  hole  through  an  aluminum  block.  The 
coil  for  the  test,  Cl,  is  shown  roughly  to  scale. 

and  then  compute  the  change  in  the  probe  impedance,  due  to  the  flaw,  as 

A Zaaw  =  Ein  ■  J f  +  Ein  •  Jint V  +  Ein  •  Jint^  (9) 

where  the  dot  products  are  the  usual  expressions  for  impedances  [1].  An  attractive 
feature  of  this  system  is  that  the  flaw  and  background  may  be  gridded  separately,  in  a 
manner  appropriate  to  their  size  and  characteristic. 

NOTCH  AT  A  BOLT  HOLE  BENCHMARK  PROBLEMS 

Benchmark  Problem  1  is  shown  in  Figure  2,  and  Benchmark  Problem  2  is  shown 
in  Figure  3.  The  two  coils  that  are  used  in  the  benchmark  problems  are  described 
in  Table  1.  Cl  will  be  used  to  collect  data  for  benchmarkl  and  HI  for  benchmark2. 
The  values  of  Cp,  L0,  and  Rq  are  determined  by  fitting  the  measured  data  of  the  coils 
in  air  with  the  usual  circuit  model  of  a  one-port  probe,  which  consists  of  a  parasitic 
admittance  element,  Yp,  in  shunt  with  a  series  impedance  i?0  +  juL0. 

MODEL-BASED  INVERSION  WITH  NLSE 

After  removing  the  shunt  parasitic  admittance  element,  Yp  (which  includes  the 
capacitance,  Cp ),  of  each  coil,  and  then  removing  the  freespace  impedance,  Rq  +  ju>L0, 
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FIGURE  3.  Benchmark  Problem  2:  Notch  in  a  bolt  hole  through  a  4mm  thick  aluminum 
plate.  The  coil  for  the  test,  HI,  is  shown  roughly  to  scale. 

TABLE  1.  Data  for  Coils  Cl  and  HI:  Dimensions  in  mm. 


Coil 

Outer  Radius 

Inner  Radius 

Height 

No.  Turns 

C'p(pF) 

L0(  H) 

Ro(to) 

Cl 

4 

1.58 

1.042 

305 

13.85 

463  x  10'6 

19.08 

HI 

7.38 

2.51 

4.99 

4000 

72.8 

0.10052 

688.4 

we  are  left  with  the  ’air-balanced’  change  in  impedance,  5Z,  when  the  coil  is  placed  over 
the  workpiece.  This  is  the  input  to  NLSE  for  inverting  each  of  the  datasets  known  as 
Ml,  M2,  M3,  and  M4  for  each  of  the  two  benchmark  tests.  (The  datasets  correspond  to 
four  different  locations  on  the  unflawed  workpiece.)  The  purpose  of  the  inversions  is  to 
determine  the  host  conductivity  and  coil  liftoff  for  each  of  the  benchmark  tests.  Even 
though  each  dataset  contained  values  of  5Z  at  32  frequencies  logarithmically  spaced 
between  500Hz  and  500KHz,  we  used  only  five  frequencies  for  the  inversions.  For 
benchmark  1,  we  used  the  five  highest  frequencies,  205.06kHz,  256.24kHz,  320.2kHz, 
400.13kHz,  and  500.0kHz,  and  for  benchmark  2  the  five  lowest  frequencies,  500Hz, 
624.8Hz,  780.76Hz,  975.65Hz,  and  1219.2Hz.  The  reason  for  this  separation  is  due  to 
the  fact  that  the  coil  for  benchmark  1,  Cl,  has  a  very  small  inductance,  as  is  seen  from 
Table  1,  whereas  that  for  HI  is  quite  large.  Therefore,  we  get  more  reliable  data  at 
high  frequencies  for  Cl  and  at  low  frequencies  for  HI.  The  inversion  results  are  shown 
in  Table  2  for  benchmark  1,  and  in  Table  3  for  benchmark  2. 

When  we  use  the  coil  data  of  Table  1  and  the  results  of  data  set  Ml  of  benchmark 
1  in  Table  2  to  compute  the  forward  model,  and  then  plot  the  results  over  the  entire 
frequency  range  of  500Hz  to  500kHz,  we  obtain  Figure  4.  (The  data  are  normalized 
with  respect  to  frequency.)  That  there  is  a  breakdown  in  the  measured  resistance  data 
at  low  frequencies  ( dR  goes  negative!)  was  anticipated  because  of  the  considerable 
uncertainty  of  Yp  for  coil  Cl  at  low  frequencies.  The  corresponding  results  for  Ml  of 
benchmark  2  in  Table  3  are  shown  in  Figure  5.  The  results  for  the  other  data  sets  in 
Tables  2  and  3  lie  on  top  of  the  results  shown  in  Figures  4  and  5.  The  high-frequency 
breakdown  in  the  resistance  and  reactance  in  Figure  5  was  anticipated  because  of  the 
uncertainty  in  Yp  for  HI  at  high  frequencies.  Note,  however,  that  the  normalized  model 
reactance,  dXj  f ,  due  to  NLSE  in  Figure  5  ’saturates’  at  high  frequencies,  as  expected 
by  classical  coupled-circuit  theory  and  Forster  plots.  The  asymptote  is  a  measure  of 
the  coupling  of  the  coil  to  the  host.  The  asymptotic  value  of  dR/ f  =  0  is  reached  at  a 
higher  frequency  in  Figure  5,  as  is  also  predicted  by  coupled-circuit  theory. 
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TABLE  2.  NLSE  inversion  results  for  benchmark  1  datasets. 


Data  Set 

$ 

cr  ( S/m)  /sensitivity 

LO  (mm)  /  sensitivity 

Ml 

4.9166 

1.8156(7)/2.73(-2) 

0. 1336/1. 262(-2) 

M2 

4.8827 

1.8151(7)/2.74(-2) 

0.1364/1. 255(-2) 

M3 

4.6704 

1.8178(7)/2.63(-2) 

0. 1420/1. 205(-2) 

M4 

4.8004 

1.8175(7)/2.69(-2) 

0. 1389/1. 236(-2) 

Notch  at  a  Bolt  Hole  Benchmark  Problem  1  Notch  at  a  Bolt  Hole  Benchmark  Problem  1 
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FIGURE  4.  Comparing  the  normalized  model  results  and  data  sets  for  Ml  of  benchmark  1 
in  Table  2  over  the  entire  frequency  range. 

TABLE  3.  NLSE  inversion  results  for  benchmark  2  datasets. 


Data  Set 

4> 

cr(S/ m) /sensitivity 

LO  (mm)  / sensitivity 

Ml 

0.6385 

1.8340(7)  /4.03(-3) 

0.3480/1. 542(-2) 

M2 

0.8585 

1.8420(7) /5.46(-3) 

0. 3560/2. 075(-2) 

M3 

0.7850 

1.8390(7)/4.97(-3) 

0.3550/1. 899(-2) 

M4 

0.3887 

1.8227(7)/2.41(-3) 

0. 3254/0. 933(-2) 

Notch  at  a  Bolt  Hole  Benchmark  Prob.  2 


Notch  at  a  Bolt  Hole  Benchmark  Prob.  2 


FIGURE  5.  Comparing  the  normalized  model  results  and  data  sets  for  Ml  of  benchmark  2 
in  Table  3  over  the  entire  frequency  range. 

Approved  for  public  release;  distribution  unlimited. 
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FIGURE  6.  Comparison  of  measured  values  of  the  change  in  impedance  to  values  computed 
using  different  grids  for  the  bolt  hole  in  Benchmark  Problem  2.  The  scan  is  along  Y  =  0. 

FINAL  RESULTS  OF  THE  BENCHMARK  TESTS 

The  final  results  of  our  validation  studies  for  the  multi-scale  algorithm  are  plot¬ 
ted  in  Figure  6  for  Benchmark  Problem  2  and  Figure  7  for  Benchmark  Problem  1. 
Numerical  experiments  show  that  accurate  results  for  Benchmark  Problem  2  require  a 
background  grid  that  is  much  finer  near  the  notch.  We  can  obtain  results  that  effec¬ 
tively  use  a  graduated  grid  which  becomes  progressively  hirer  in  regions  closer  to  the 
notch.  It  is  clear,  however,  from  varying  the  number  of  notch  cells  that  the  impedances 
we  are  trying  to  compute  have  not  converged  with  respect  to  number  of  cells  along  the 
x,  y,  or  z  directions  for  a  grid  of  32  x  2  x  16  cells. 

It  is  also  clear  from  varying  the  number  of  background  (bolt-hole)  cells  that  the 
impedances  have  not  converged  with  respect  to  number  of  cells  along  the  x,  y,  or  z 
directions  for  a  grid  of  16  x  16  x  4  cells. 

We  can  improve  our  impedance  values  by  combining  them  with  values  obtained 
from  an  8  x  8  x  2  cell  graduated  grid  for  the  bolt  hole  and  a  16  x  2  x  8  cell  notch  grid. 
This  allows  us  to  perform  a  linear  extrapolation  to  zero  cell  dimensions.  That  is  the 
result  labeled  ’Extrapol’d’  in  Figure  6. 

As  with  Benchmark  Problem  2,  accurate  results  for  Benchmark  Problem  1  re¬ 
quire  a  background  grid  that  is  much  hirer  near  the  notch,  and  so  we  again  apply  the 
graduated-grid  scheme  that  was  used  for  that  problem.  We  also  use  the  extrapolation 
procedure  that  was  used  in  Problem  2,  with  the  result  labeled  ’Extrapol’d’  in  Figure  7. 

The  response  of  the  small  notch  is  clearly  evident  in  both  benchmark  problems, 
and  its  shape  agrees  with  the  measured  data,  which  were  the  goals  of  the  algorithm,  but 
the  absolute  values  differ  by  12%  to  39%,  with  the  smaller  errors  occuring  in  Benchmark 
Problem  2.  The  reason  for  the  discrepancies  in  scale  is  due  to  the  size  of  the  problems. 
Benchmark  Problem  1  contains  a  hole  that  is  abnormally  large  compared  to  the  notch, 
and  at  the  excitation  frequency  of  5kHz,  this  hole  requires  a  huge  number  of  cells. 
Benchmark  Problem  2  is  somewhat  more  realistic  in  its  size,  and  the  extrapolated 
results  are  closer  in  scale  to  the  measured  data.  We  have  ideas  for  improving  the 
convergence  of  the  algorithm  for  large  problems,  and  they  will  be  explored  in  future 
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FIGURE  7.  Comparison  of  measured  values  of  the  change  in  impedance  to  values  computed 
using  different  grids  for  the  bolt  hole  in  Benchmark  Problem  1.  The  scan  is  along  Y  =  0. 


research. 
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